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Abstract. This research describes a three dimensional quantum cellular automa- 
ton (QCA) which can simulate all other 3D QCA. This intrinsically universal QCA 
belongs to the simplest subclass of QCA: Partitioned QCA (PQCA). PQCA are 
QCA of a particular form, where incoming information is scattered by a fixed uni- 
tary U before being redistributed and rescattered. Our construction is minimal 
amongst PQCA, having block size 2x2x2 and cell dimension 2. Signals, wires 
and gates emerge in an elegant fashion. 



1. Introduction 

(Quantum) Cellular Automata. Cellular automata (CA), first introduced by von 
Neumann [52j, consist of an array of identical cells, each of which may take one of 
a finite number of possible states. The entire array evolves in discrete time steps by 
iterating a function G. This global evolution G is shift-invariant (it acts the same 
everywhere) and causal (information cannot be transmitted faster than some fixed 
number of cells per time step). Because this is a physics-like model of computation 
[27j, Feynman [22j, and later Margolus [28], suggested early in the development 
of quantum computation (QC) that quantising this model was important. This 
was for two main reasons. Firstly, in QCA computation occurs without extraneous 
(unnecessary) control, hence eliminating a source of decoherence, which is problem 
for QC implementation. Secondly, they are a good framework in which to study 
the quantum simulation of a quantum system, which is predicted to be a major 
application of QC. From a Computer Science perspective there are other reasons 
to study QCA, such as studying space-sensitive problems in computer science {e.g. 
'machine self-reproduction'. Tiring Squad Synchronisation', ...) in the quantum 
setting, or having more complete and formal models of physics ability to compute. 
There is also the theoretical physics perspective, where CA are used as toy models 
of quantum space-time [25] . The first approach to defining QCA [2[ [191153] was later 
superseded by a more axiomatic approach |9l [IHl [J4J together with more operational 
approaches [I5l [36l SH [42l [49l . 
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Intrinsic universality. Probably the most well known CA is Conway's 'Game of Life', 
a two-dimensional CA which has been shown to be universal for computation, in 
the sense that any Turing Machine (TM) can be encoded within its initial state 
and then executed by evolution of the CA. As TM are often regarded as a robust 
definition of 'what an algorithm is' in classical computer science, this provides a 
key result in CA research. However, more can be achieved using CA than just 
running any algorithm. They run distributed algorithms in a distributed manner, 
model phenomena together with their spatial structure, and allow the use of the 
spatial parallelism inherent to the model. These features, modelled by CA and not 
by TM, are all worthy of investigation, and so the concept of universality should 
be revisited in this context to account for space. This is achieved by returning to 
the original meaning of the word universality [H [HI [16] , namely the ability for one 
instance of a computational model to be able to simulate other instances of the 
same computational model. Intrinsic universality formalises the ability of a CA to 
simulate another in a space-preserving manner [29l |39l [47] , and has been extended 
to the quantum setting [4[ |6l [7]. In previous work [7] it was shown that, when the 
dimension of space is greater than one, the problem of intrinsic universality reduces 
to the ability to code for signals, wires and a universal set of quantum gates, as 
expected from the classical C A case [39j . 

Related work. There are several related results in the current CA literature. For 
example, several works [271 EH |35j provide computation universal Reversible Parti- 
tioned CA constructions, and deals with their ability to simulate any CA in the 
one-dimensional case. The problem of minimal intrinsically universal CA has been 
addressed ^40], and for Reversible CA (RCA) the issue of intrinsic universality has 
been tackled flTl dU [23] . The difficulty is in having an n-dimensional RCA simulate 
all other n-dimensional RCA and not, say, the (n— l)-dimensional RCA, otherwise a 
history- keeping dimension could be used, as shown by Toffoli [48]. Concerning QCA, 
Watrous proved that QCA are universal in the sense of QTM. Shepherd, Franz 
and Werner [45] defined a class of QCA where the scattering unitary Ui changes at 
each step i (CCQCA). Universality in the circuit-sense has been achieved by Van 
Dam ^], Cirac and VoUbrecht [51j, Nagaj and Wocjan [36j and Raussendorf [42j. 
In the bounded-size configurations case, circuit universality coincides with intrinsic 
universality, as noted by Van Dam |49j. QCA intrinsic universality in the one- 
dimensional case has been resolved [5], and also in the general n > 1-dimensional 
case [7J . Both results rely on recent work [6j , where it was shown that a simple sub- 
class of QCA, namely Partitioned QCA (PQCA), are intrinsically universal. This 
enables the focus here to be on this simple, natural class of QCA, as previously 
proposed [13 [SI |36l [El [H [IS [53] . PQCA are QCA of a particular form, where 
incoming information is scattered by a fixed unitary U before being redistributed 
and rescattered. Hence the problem of finding an intrinsically universal PQCA is 
reduced to finding some scattering unitary U that supports the implementation of 
signals, wires and a universal set of quantum gates. 

Game of Life. As the more general case of finding a general n-dimensional intrin- 
sically universal QCA has been resolved [i6], one could question the necessity of 
focusing on the particular 3D case of this problem. Our answer to this question is 
threefold. 
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• In the study of models of computation, the definition of the model is often 
followed by the search for a universal instance, which is then followed by a 
search for a minimal universal instance. Reaching this final step is generally 
regarded as a sign of maturity of the model. It shows that the entire model 
can be reduced to a particular simplest instance. Moreover, sometimes a 
particular instance stands out, either for its elegance and simplicity, or the 
richness of its behaviour. This was clearly the case with the Game of Life 
within the realm of Classical CA, and it is echoing this impression that we 
have named this QCA the 'Quantum Game of Life'. This is not because of 
a resemblance in the construction of the local rule, which clearly owes more 
to the Bilhard BaU Model CA [23]. 

• Dimension three is particularly interesting, not only because of the relevance 
to physics, but because intrinsic universality QCA constructions lend them- 
selves to a striking simplification. In dimension one the simplest known 
intrinsically universal QCA has cell dimension 36 [5j. In dimension two it 
has cell dimension 4 [7J , and seems reducible to 3 via a costly and inelegant 
variation of the scheme [7j, but probably no further. The solution we present 
here is in dimension three and has cell dimension 2; it is therefore minimal. 

• Animations of the QCA operating, plus an implementation of this QCA, are 
available on the companion website [8j. Although it lacks a friendly way of 
preparing the initial configuration, it turns out to be quite interesting and 
entertaining already, as is the Game of Life. We plan to implement the QCA 
in a more graphical, interactive, computationally efficient fashion. The end 
result could be a pedagogical tool, both for discussing quantum theory and 
in order to start exploring the dynamical behaviour of QCA, such as in [37] . 
perhaps using methods inspired by similar work on Probabilistic CA [43j. 

2. The Universal QCA 

Definitions. Configurations hold the basic states of an entire array of cells, and hence 
denote the possible basic states of the entire QCA: 

Definition 1 (Finite configurations). A (finite) configuration c over S is a function 
c : — > S, with (^1,^2, ^3) ' — > ^(^1,^2, ^3) = ^1,^2,^35 ^^^^ that there exists a 
(possibly empty) finite set / satisfying (^1,^2, ^3) ^ I ^ ^hMM ^ where g is a 
distinguished quiescent state of S. The set of all finite configurations over S will be 
denoted 6^. 

As this work relates to QCA rather than classical CA, the global state can be 
a superposition of these configurations. To construct the separable Hilbert space of 
superpositions of configurations the set of configurations must be countable. This 
is why finite, unbounded, configurations are considered; the quiescent state of a CA 
is analogous to the blank symbol of a TM tape. 

Definition 2 (Superpositions of configurations). Let IKgs be the Hilbert space of 

configurations. Each finite configuration c is associated with a unit vector |c), such 
that the family (|c))^^gs is an orthonormal basis of IKgs. A superposition of config- 
urations is then a unit vector in IKps . 
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Figure 1: Signals propagate diagonally across the partition. 

Detailed explanations of these definitions, as well as axiomatic definitions of 
QCA, are available [9l [IHl E]. Building upon these works, we have shown [61 17j 
that Partitioned QCA (PQCA) are intrinsically universal. Since they are the most 
canonical description of QCA, and since the aim of this paper is to construct the 
most canonical, and yet universal, QCA, we will assume that all QCA are PQCA 
throughout this work. 

Definition 3 (Partitioned QCA). A partitioned three-dimensional quantum cellular 
automaton (PQCA) is defined by a scattering unitary^ a unitary operator, U such 
that U : IKf — > 3^1^^, and U \ qq . . . qq) = \qq . . . qq)^ i.e. that takes a cube of 2^ 
cells into a cube of 2^ cells and preserves quiescence. Consider G = {^21? 
operator over "K. The induced global evolution is G at odd time steps, and aG at 
even time steps, where a is a translation by one unit in all directions. 

Here we provide a specific instance of a [/-defined PQCA which is capable of 
intrinsically simulating any ^/-defined PQCA, for any V . In order to describe such 
a [/-defined PQCA in detail, two things are required: the dimensionality of the cells 
(including the meaning attached to each of the states they may take); and the way 
in which the scattering unitary U acts upon these cells. 

2.1. Signals and wires 

Classical CA studies often refer to 'signals' without an explicit definition. In 
this context, a signal refers to the state of a cell which may move to a neighbouring 
cell consistently, from one step to another, by the evolution of the CA. Therefore a 
signal would appear as a line in the space-time diagram of the CA. 

The [/-defined PQCA constructed in this paper is minimal, in the sense that it is 
a 3D PQCA with the smallest non-trivial block-size and cell-dimension. Hence, each 
cell has only two possible states, either empty or occupied; it is a binary automaton, 
with scattering unitary U acting on 2 x 2 x 2 cell neighbourhoods (a cube). The 
description of the scattering unitary U that defines the behaviour of the automaton 
will now be provided in a case by case manner. 

2.1.1. Signals. Signals in our scheme will be encoded as isolated occupied cells. 
If a neighbourhood trivially contains no signals, then no change is made. If a 
neighbourhood contains only one signal, i.e. only one cell is occupied, and the other 
seven are empty, then the signal is shifted by one in every direction, yielding a 3D 
diagonal propagation. This rule is given in Fig. [T] These two rules, together with 
their rotations, define the action of the scattering unitary U upon the subspace with 
zero or one occupied cells. Notice that U can easily be seen to be unitary upon this 
subspace, as it performs a permutation of the basis states. 
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Figure 2: Two occupied cells which share one face form a static barrier. 




Figure 3: Two vertical barriers placed in adjacent partitions horizontally form a 
2x2x1 barrier (left). At the next time step, the repartitioning is still 
stable, with two horizontal barriers in adjacent vertical partitions formed 
from those same cells (right). 

$ 

Figure 4: Four occupied cells occupying a face of the partition neighbourhood form 
a wall. 

2.1.2. Barriers. Two or more occupied cells that share a common face form a barrier, 
which is intended to be a stationary pattern. The rule is given in Fig. [2| Such 
a barrier is stable for one partition, but in the next partition it appears as two 
signals moving away from each other, and hence it will scatter. This can be avoided 
by extending beyond the current partition, using at least four occupied cells. For 
example, a 2 x 2 x 1 block of cells which overlaps a partition and forms a barrier in 
each partition, as shown in Fig. [3| is stable. This rule, together with its rotations, 
defines the action of the scattering unitary U upon the subspace of two adjacent 
occupied cells. Upon this subspace the scattering unitary U acts like the identity, 
which is unitary. 

2.1.3. Walls. When two barriers, made of four occupied cells, form a square that fits 
within a partition, as in Fig.[4| we say that they form a wall. A wall is able to redirect 
a fifth occupied cell, that shares only one face with the wall, which represents the 
incoming signal. Again a wall in one partition appears as signals moving away from 
each other in the next partition, but again they can be made stable following Fig. [3| 
by extending them across partitions. The behaviour required is for an incoming 
signal to "bounce" off the wall, and the rule producing this behaviour is given in 
Fig. [5| This rule suflBces to define signal redirection, or "rewiring". These rules, 
together with their rotations, define the action of the scattering unitary U upon 
the subspace of walls and walls with a signal. Notice that U is unitary upon this 
subspace, as it performs a permutation of the basis states. 
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Figure 5: A signal which interacts with a wall in a neighbourhood "bounces" off 
the wall, causing it to change direction along one axis. Compare with the 
unhindered signal propagation shown in Fig. [TJ 




Figure 6: Qubits are implemented as parallel tracks of signals. No qubit is modelled 
as no signal (/e/i), |0) is modelled as a signal on the lower track (centre)^ 
while a signal on the upper track models a qubit in state |1) (right). 



4; 





















Figure 7: The rules implementing the Hadamard operation. The first gives |0) 
^ |0) + ^ |1) while the second gives |1) ^ 73 |0) - ^73 

2.2. Qubits and quantum gates 

To allow a universal set of gates to be implemented by the PQCA, certain com- 
binations of signals and barriers need to be given special importance. To implement 
qubits, pairs of signals will be used. Indeed in our scheme a qubit is formed by 
two parallel "tracks" of signals, as shown in Fig. |6| A signal on the bottom track 
indicates a |0) state, whereas the top track indicates a |1) state. Qubits in super- 
positions are modelled by appropriate superpositions of the configurations, as in 
Def . § We can now define rules which implement quantum gates on those qubit^ 



2.2.1. Hadamard. The Hadamard gate is a requirement for universal quantum com- 
putation, where i7 : |0) 1-^ |0) + |1) , |1) |0) — In order to achieve 
this typically quantum behaviour, a special meaning is attached to the interaction of 



a signal with a single two-cell barrier (as defined in section 2.1.2). The corresponding 
rule is given by Fig. [7| 

■'^In Quantum Mechanics, a qubit can be encoded into any 2 degrees of freedom of a system. 
Some experiments in quantum optics, for example, use the spin degree of freedom, whereas others 
use a spatial degree of freedom. This model is analogous to the latter case. See the conclusion 
section for a discussion of relevance to physical systems. 
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Figure 8: A stable Hadamard configuration, formed by two barriers extended across 
partitions (left). The Hadamard configuration can be extended with two 
cubed "waUs". The first cube, on the top left, redirects the top track |1) 
signal, incoming from the bottom left and travelling to the far top right, 
into the barrier. The second cube, bottom right, redirects the output 
bottom |0) track so that it is again moving parallel to the |1) track, keeping 
the qubit representation consistent (right). 

As all the rules we provide are rotation invariant, the unambiguity of the Hadamard 
rule may not be clear. To be unambiguous, the |0) and |1) cases need be to distin- 
guished from each other, and the signal distinguishable from the barrier. This can 
be done as follows: if the cube can be rotated such that all three signals are on a 
single face, forming an "L" shape, then the input to the Hadamard is |0), and the 
signal is at the top of the L. If instead a lone occupied cell is on the top left of the 
furthest face, with a two-cell barrier along the bottom of the closest face, forming 
a dislocated L shape, then the input to the Hadamard is |1) and the signal is the 
isolated, top left cell. 

On its own, a barrier is not enough to implement the Hadamard gate, because 
it does not respect the parallel track convention of qubit s defined earlier. Rather, 
inputs and outputs are orthogonal. To rectify this, the barrier can be combined 
with two walls, which redirect the signals appropriately, as shown in Fig. [8| Further 
rewiring can be added so the qubit can arrive and leave in any direction, but these 
are not easily represented in 2D projections. 

These two rules, together with their rotations, define the action of the scattering 
unitary U upon the subspace of two-cell barriers plus a signal. Notice that U is 
unitary upon this subspace, as it is a composition of a Hadamard gate upon the 
input states, together with the permutation of the basis states that perform the 
movement described in Fig. [7| 

2.2.2. Controlled rotations. A two qubit controlled gate is another requirement for 
universal quantum computation, and in this case we choose the controlled-R(|)gate, 

where cR(7r/4): |11) |11), and is the identity otherwise. To encode a two 

qubit controlled gate, signal collisions, where two moving signals cross each other 
diagonally, will be interpreted as adding a global phase of e^. This rule is given by 
Fig-i 

This allows a controlled-R(|) operation to be defined, by redirecting the |1) 
(true) signal track so that the tracks will cross each other in this way, and then 
recreating the qubits by rejoining the |1) tracks back with the appropriate |0) tracks, 
using walls to redirect and delay signals as in the Hadamard case. Again, this 3D 
configuration cannot easily be presented in 2D. 
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Figure 9: When two signals cross each other diagonaUy, a complex phase is added 
to the configuration state. 



This same rule can also be used to implement a single qubit phase change op- 
eration by causing a control signal to loop such that it will intersect with the true 
track of the qubit, thus adding the global phase to the configuration. This needs 
to be correctly timed so that the tile implementing the single qubit rotation can be 
iterated, following the same reasoning as in [7J. 

This rule defines, together with its rotations, the action of the scattering unitary 
U upon the subspace of two non- adjacent signals on face. U can again be seen to be 
unitary upon this subspace, as it is a permutation of the basis states, with a phase. 



2.3. Universality, erosion and dynamics 

We have seen that the scattering unitary [/ is a unitary on the subspaces upon 
which it has been defined. Additionally, following arguments given previously [7], it 
may be concluded that the [/-defined PQCA given here is intrinsically universal: 

• In space dimension greater than one, it suflBces to implement signals, wires 
and a universal set of gates to achieve intrinsic universality; 

• The controlled-R(l) and the Hadamard gate are universal for Quantum 
Computation. 

Indeed the standard set of cNOT, H, R(|) can be recovered as follows: 

CNOT = (I ® //)(cR(7r/4))^(][ ® H) \^) 

where cR(|)^ denotes four applications of the controlled-R(|) gate, giving the 
controUed-PHASE operation. 

We could let the action of the scattering unitary U be the identity in all other 
subspaces. However, some patterns would then be unmovable and indestructible, 
such as the wall of Fig. [4] This would arguably make the dynamics of the QCA 
presented here a little plain. Moreover, since walls could be built through simulta- 
neous signal collision, it would be natural to be able to dismantle them. If we are 
interested in the dynamics, diflFerent ways of completing the definition of U upon 
the remaining subspaces should be considered. A natural way is to consider that 
all of the other cases are made of non-interacting signals, where the occupied cells 
propagate past each other diagonally. This allows, for instance, for the demolition 
of walls and the creation of new stable patterns. 

Animated examples and an implementation of this cellular automaton can be 
found online [8]. We hope to explore its dynamics in future simulations. 



3. Comparison with Two-Dimensional Rules & Tiles 

In previous work we provided a generic construction of an n-dimensional in- 
trinsically universal QCA pj. In this section we will compare the two dimensional 
instance given previously with the scheme presented here. 
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Figure 10: The 2D 'controlled-R(|) gate' tile, with a signal interaction at the high- 
lighted cell. 



Objects. In [7] the cell-dimension is 4. A cell can be either empty, a wall, or a 
binary signal (0 and 1). In a sense the same objects are recovered in the scheme 
presented here, but as the cell dimension is 2 these are conveyed as patterns of cells. 
Walls are now coded as 2^ squares as in Fig. [ij Signals no longer have a bit attached 
to them; instead the location of the signal is used to encode this bit of information. 

Rules. Once the above observations have been made; it can be seen that the 
two dimensional, cell-dimension four scheme previously described [7J is almost a 
projection of the scheme presented here. Indeed: 



- the 2D signal travelling rule of [7j; 
of the corresponding rule in Fig. [T] 

- the 2D signal bouncing rule; 
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the corresponding rule in Fig. 
- the 2D control-phase rule; 
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is a width-depth projection of the corresponding rule in Fig. |9| 
- only the 2D Hadamard rule of |2J differs slightly: 
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Tiles. Previously, not only were the above rules provided, but a set of fixed-sized 
(16 X 14) tiles, taking a fixed number of time steps (24) to be travelled through, and 
each accomplishing one of the universal quantum gates upon the incoming signals 
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were also given [7J. The example of the controlled-R(|) tile is reproduced in Fig. 
for concreteness. The difference between such tiles and the corresponding rule (the 
controlled-R(l) rule in this case) is that tiles need only to be placed next to one 
another so as to make up a circuit, with all signal wiring and synchronisations being 
accounted for. It would be cumbersome to describe the equivalent of all these tiles in 
the scheme presented here, particularly due to their three dimensional nature. Some 
are shown on the companion website [8j . The reader may convince himself of their 
feasibility through the following arguments: since those tiles exist in two-dimensions, 
and since the two-dimensional rules are projections of the three-dimensional rules. 
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there must exist some three-dimensional tiles for which the two-dimensional ones 
are, in some sense, projections. 

4. Conclusion 

This paper presents a minimal 3D PQCA which is capable of simulating all 
other PQCA, preserving the topology of the simulated PQCA. This means that the 
initial configuration and the forward evolution of any PQCA can be encoded within 
the initial configuration of this PQCA, with each simulated cell encoded as a group 
of adjacent cells in the PQCA, i.e. intrinsic simulation. The main, formal result of 
this work can therefore be stated as: 

Claim 1. There exists an 3D [/-defined PQCA, with block size 2 and cell dimension 
2, which is an intrinsically universal PQCA. Let be a 3-dimensional ^-defined 
PQCA such that V can be expressed as a quantum circuit C made of gates from 
the set Hadamard, Cnot, and R(|). Then G is able to intrinsically simulate H. 

Any finite-dimensional unitary V can always be approximated by a circuit CiV) 
with an arbitrary small error e = max|^) \ \V — C 1 1. Assuming instead that G 
simulates the (7(l/)-defined PQCA, for a region of s cells over a period t, the error 
with respect to the ^/-defined PQCA will be bounded by ste. This is due to the 
general statement that errors in quantum circuits increase, at most proportionally 
with time and space [38j. Combined with the fact that PQCA are universal [TJ [6], 
this means that G is intrinsically universal, up to this unavoidable approximation. 

Discussion. This PQCA is definitely minimal amongst the 3D PQCA. Reducing 
the block size or the cell dimension necessarily results in a trivial PQCA. Moreover, 
PQCA are the simplest and most natural class of QCA [3j. Nevertheless, it is not so 
clear that this PQCA is minimal amongst all intrinsically universal 3D QCA. Indeed, 
PQCA-cells are really QCA-subcells, and PQCA-blocks need to be composed with 
them shifted in order to yield the QCA neighbourhood. Based on these observations 
our QCA has QCA-cell dimension 2^ and the 3D radius-| neighbourhood. Whether 
there are some further simplifications in this more general setting is an open question. 

An important source of inspiration, along with the original Game of Life, was 
the the simplest known intrinsically universal classical Partitioned CA [28], which 
has cell dimension 2. Called the BBM CA, it was itself directly inspired by the Bil- 
liard Ball Model |23j. Our scheme also features signals (billiard balls, or particles) 
bouncing and colliding. This analogy with physics is a desirable feature; Feynman's 
sole purpose for the invention of QC [21 J was that quantum computers would be able 
to simulate quantum systems efficiently, and hence predict their behaviour. It is still 
thought that quantum simulation will be one of most important uses of quantum 
computers for society, with expected and potentially unexpected impact in quan- 
tum chemistry, biochemistry or nanotechnologies {e.g. in terms of the synthesis of 
specific purpose molecules). This was also one of the reasons why Feynman immedi- 
ately promoted the QCA model of QC [22]; they are a natural mathematical setting 
in which to encode the continuous-time and space, sometimes complex, behaviour 
of elementary physical elements (such as wave equations of particles, possibly in- 
teracting, under some field, for example) into the more discrete-data discrete-time 
framework of quantum computation. This issue of encoding is non-trivial, but has 
largely been solved for the one-dimensional case, by using QCA as a mathematical 
framework in which to model the system of interest. For example, several works 
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[lailSlIHEQlEHlEOlEIllaS^ simulate quantum field theoretical equa- 

tions in the continuous limit of a QCA dynamics. These results do not extend to 
more-than-one dimensions due to the problem of isotropy. Whilst in one-dimension 
the division of space into a line of points is an acceptable model, because it does 
not privilege the left over the right direction, in two-dimensions the grid seems 
to inevitably favour the cardinal directions (N, E, S, W) over the diagonal ones. 
That this construction is similar to the Billiard Ball Model CA suggests a Quantum 
Billiard-Ball Model could be defined, moving away from the underlying grid, which 
would remedy this problem. This is planned this for future work. 

After this work was completed, a two-dimensional, continuous-time, quantum 
version of the Game of Life was reported (see arXiv:1010.4666 ) . 
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